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ABSTRACT 

Neutron-induced reaction rates, including fission and neutron capture, are calculated in the temperature range 10 s < T(K) < 10'° 
within the framework of the statistical model for targets with the atomic number 84 < Z < 118 (from Po to Uuo) from the neutron 
to the proton drip-line. Four sets of rates have been calculated, utilizing - where possible - consistent nuclear data for neutron separa- 
tion energies and fission barriers from Thomas-Fermi (TF), Extended Thomas-Fermi plus Strutinsky Integral (ETFSI), Finite-Range 
Droplet Model (FRDM) and Hartree-Fock-Bogolyubov (HFB) predictions. Tables of calculated values as well as analytic seven pa- 
rameter fits in the standard REACLIB format are supplied . We also discuss the sensitivity of the rates to the input, aiming at a better 
understanding of the variations introduced by the nuclear input. 
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1. Introduction 

Investigations of nucleosynthesis processes make use of reaction 
networks including thousands of nuclei and tens of thousands of 
reactions. Most of these reactions occur far from stability and 
thus cannot yet be directly studied in the laboratory. In addition 
most of the nuclear properties including reaction rates, which 
are also required for the calculation of cross sections and as- 
trophysical reaction rates, are not experimentally known either. 
Therefore, predictions based on theoretical models are neces- 
sary. While close to stability partial experimental information is 
available, relying fully on theoretical information leads to rela- 
tively large variations in computed cross sections far from sta- 
bility. This is especially true for the region of fissionable nuclei, 
which is the focus of the present investigation. 

In the past, a series of efforts were applied to calculate 
neutron-capture rates for r-proc ess nucleosynth e sis and other 
astrophysical applications ( e.g.. | Arnouldl fl972t iHolmesetaEI 
197fj[ IWooslev et ail 1 19781 ISargoodl Il982t IThielemann et all 
19871 ICowan et alj 11991b IRauscher & Thielemann! 12000; 



Aika wa et al.ll2005 : Gorielv et al.ll2008l and references therein - ). 

Fission has often been neglected in astrophysical calcu- 
lations. In early applications to astrophysical nucleosynthesis, 
usually only one mode was considered, beta-delayed fission 
(IThielemann et a l l 19831) or a phenomenological model of spon- 
taneous fission ([Go i'ielv & Clerbauxl Il999t iFrei burghaus et alj 
U999HCowanetalJll999h . However, it was shown recently that 
neutron-induced fission is more im portant than beta-delayed fis- 
sion in r-process nucleosynthesis dPanov & Thielemann! 120031 



12004 iMartfnez-Pinedo et alj 120071) . Thus, the need to pro- 
vide a compilation of neutron-induced fission ra tes is obvi- 
ous. In itial investigations have been undertaken by iPanov et"a"D 
d2005l) and lGorielv et alj d2009l) . Here we present extended cal- 
culations of neutron-induced fission rates for different predic- 
tions of masses and fission barriers. The present work also com- 
pletes existin g nuclear neutron-captur e rate sets b y extending 
the wo rks of IRauscher & Thielemann! d2000t) and IPanov et"al] 
(2005) to the region 84 < Z < 118 in order to provide the 
necessary input for nucleosynthesis stu dies under high neu- 
tron densities. As in IPanov et"aT1 (120051) . t he statistical model 
approach of Wolfenstein- Hauser-Feshbach dWolfensteinl [19511 
lHauser & Feshbach ll952l) for compound nuclear reactions was 
used, but employing more recent and complete data and predic- 
tions for masses, spins, and fission barriers. 

Nuclear mass and fission barrier predictions have a strong 
model dependence, and none of the existing models can repro- 
duce all experimentally known data. Moreover, the fission pro- 
cess itself is complicated, and extended calculations for neutron- 
induced fission across the nuclear chart have to be done care- 
fully. Here, we aim to provide rates for studying the endpoint 
of the r-process and the possible production of super-heavy el- 
ements. By comparing rates obtained with different choices of 
mass and fission barrier predictions we attempt to give a measure 
of the involved variations. Astrophysical models, providing the 
nucleosynthesis conditions, bear large variations in themselves. 
This is especially true for the r-process, for which the astrophys- 
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Tables 3-18 with these data are only available in electronic form at 
the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 
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ical site is still unknown despite decades of study. For a realistic 
and exhaustive exploration of synthesis conditions, simulations 
do not only have to vary astrophysical parameters, but also have 
to include a variation range of involved reaction rates given by 
different mass and fission barrier models. 

Our paper is structured as follows. In Sect. |2]we briefly de- 
scribe the statistical model used in the calculations as well as 
the nuclear input data and give a comparison of cross sections 
or rates for a number of experimentally known nuclei with ex- 
isting experimental information and other theoretical models. 
These methods are then appl ied to supplement the rate sets of 
iRauscher & Th ielemann (2000) of (n,y)-rates for chemical ele- 
ments with Z > 83 and predict neutron-induced fission cross sec- 
tions and rates (where available in comparison to experiments). 
Section [3] presents these results and shows the sensitivity with 
respect to mass models and fission barriers employed. Rate fits 
for utilization in astrophysical calculations are discussed in Sect. 
|4] In Sect.[5]we give a brief discussion and some examples of the 
mass distribution of fission fragments, which will be provided in 
an extended way in a forthcoming paper. The final Sect. [6] con- 
tains conclusions and a summary. The explanation of the tables 
and their structure are given in Appendix A. The complete ta- 
bles of reaction rates and their fits are found at CDS in electronic 
form. 



2. The statistical model and nuclear data input 

As in previous approaches (e.g.. IThielemann et ail Il987b 
ICowan et alj|199lt IRauscher & Th ielemann COOOl) we have ap- 
plied the statistical Wolfenstein-Hauser-Fe shbach formalism 
<IWorfensteinlll95lb lHauser & Feshbachlll952h for the calcula- 
tion of neutron-induced cross sections and reaction rates. In ad- 
dition to (» , y)-reactions the f i ssion channe l was also included a s 
outlined in IThielemann et al.1 (1 19831 fl 9891) . ICowan et all (I 199 ll) . 
iPanov et alj (120051) . The statistical model is applicable for as- 
trophysical rate calculations as long as there is a sufficiently 
high density of excited states in the compound nucleus at the 
relevant bombarding energy, which is the case for most heavy 
nuclei. However, at shell closures and with decreasing neutron 
separation energies, level densities at the astrophysically rele- 
vant compound formation energy in neutron-induced reactions 
be come too small for th e application of the model, as shown 
bv IRauscher et al.l dl997l) . In those cases, single resonances and 
contributions from the direct reaction mechanism have to b e 
taken into account (IRauscher et al.ll998blG"orieiv & Khanl2002l) . 
This underlines on the one hand that reliable mass predictions 
for the separation energies are absolutely necessary for r-process 
applications far from stability, and on the other hand that the 
prediction of spectroscopic properties for resonant and direct 
capture awaits improvement as well. Here we calculate neutron- 
induced rates still based purely on the statistical model to pro- 
vide a full set of rates for extended r-process calculations and 
the possible formation of superheavy nuclei for a variety of dif- 
ferent sets of mass models and fission barrier predictions. The 
influence of the direct reaction mechanism on the rates far from 
stability needs to be explored separately in a future study. 

We outlined the general treatment within the statistical 
model for applications, the cross section for a neutron-induced 
reaction f(n,out) ("out" standing for gamma-emission or fis- 
sion) from the target ground state f with center of mass energy 
E in and reduced mass fi in given by 



2(2/+ 1) 



In 



T^(E,J",E°,jf,7^)T out (E,f) 
TUE,J") 



(1) 



(n.out) 



nh 2 



1 



(2(x in E in ) (2jf + 1) • (2J n + 1) 



The total transmission coefficient T tol = £ vo T v a describes the 
transmission into all possible bound and unbound states v in 
all energetically accessible exit channels o (including the en- 
trance channel z). The fission transmission coefficient Tf(E, J*) 
includes the su m over all possible final s t ates and is evaluate d 
as discussed in iBiornholm & Lvnnl d!980l) . ICowan et alj d!991l) . 
Pano v et a D (12001 " and is related to the fission probability 
Pf(E, J") = Tf(E, J^/T ptjE, J K ) consider ed in the papers cited 
above. Since the work of Strutinskyl d 19671) fission has been gen- 
erally described within the framework of double-humped fission 
barriers. When making use of a double-humped fission barrier, 
the fission transmission coefficients can be calculated in the limit 
of complete damping, which averages over transmission reso- 
nances, assuming that levels in the second minimum between 
the first and second barrier E& and Eg are equally spaced. If de- 
tailed information concerning the level structure in the second 
well is missing, this method gives the best results. Pj further re- 
quires the transmission coefficients through the first and second 
barrier Ta and Tg, which are evaluated by an integral over the 
first and second barrier potential, weighted with the level density 
at the appropriate energy and corresponding deformation. The 
individual barriers can be approximated by individual (inverted 
parabola) Hill-Wheeler barrier shapes. The level densities p& 
and pb show an enhancement over the level densities at ground 
state deformation, and it is important to include proper symme- 
try classes in the calculations at corresponding saddle points. In 
the absence of detailed information, standard factors of four (for 
the axially asymmetric/mass symmetric barrier) and two (axially 
symmetric/mass a symmetric) over the grou nd state level den- 
sity were applied feiornhol m & Lvnnlll980l) . Whenever possi- 
ble, experimentall y known fission barri ers were used, taken from 
ISmirenkinl d 19931) . t he compilation of IMamdouh et all dl998l) . 
and the database of iBelgva et al.l (l2006|). The other ttansmis- 
sion c oefficients were calculated as in IRauscher & Thielemannl 
(2000), utilizing up to 19 experimen tally known excited states (i f 
available). The data were taken from F irestone & Shirlevl dl996). 
up to the first level for which the spin assignment is unknown. 
Ground state spins and parities are known for many unstable nu- 
clei. Far off stabili t y, gro und state spins and parities were taken 
from Mol ler et al.l (1 19971) when experimental values were un- 
av ailable. Above t he las t known state, the nuclear level density 
of Rauscher et al] d 1997b was used. This method is based on the 
back shifted Fermi-gas approach, where the level density param- 
eter a and back shift 5 are obtained globally from the appropriate 
mass model employed. 

In Fig. Q] we compare our predictions for neutron-induced 
fission cross sections of some U, Np a nd Pu isotopes 
with e valuated neutron data from JENDL-3.3 (fNakagaw a et al.l 
12005b ISonnera et al.l 120081) . (The accuracy of the evalu- 
ated data is usually not declared. In the region of interest, 
0.01 < < 10 MeV, the accura cy of up-to-dat e mea- 
surements for plutonium isotopes by iTovesson et al.l ([2009) 
varies from 2% to 15%, and for our plots is not bigger 
than the plot signs. The detailed experimental information 
can be found in the experimenta l nuclear reaction d ata library 
(2009)). Experimental masses and fission barriers or masses 
and fission ba rriers were employed from different mass mod- 
els: ETFSI dAboussir et all 1 1995b IMamdouh et ail 1 19981). TF 
dMvers & Swiateckill 1996b 1 19991) . HFB-14 dGorielv et al.ll2 009) 
and th e older liquid drop predictions bv iHoward & M oller 
(119801) . here shown as Panov et al. 2005). It can be seen that 
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Fig. 1. Present predictions of energy-dependent («,/) cross sections cr n f{E) for some target nuclei of U, Np and Pu calculated 
in the framework of different mass and fission barrier predictions (ETFSI, TF, HFB-14) and exp erimental data, marked Bj Xp as 
well. Experimentall y measured cross-sections were used after JENDL-3.3 (Nakagaw a et alj|2005l) . averaged by the code JAMS 
Sopp era et al] ( 20081) . displayed by a black line. All the predictions are given for a ground-state population. Our previous results 
dPanov et al.ll2005l) are shown as well. 



when using experimental fission barriers, the agreement with 
experimental («,/) cross sections is within the typical factor 
of two to three known for statistical model calculations. It 
can also be seen that different barrier predictions can lead to 
large variations, which will clearly remain for predictions far 
from stability where no experimental information is available. 
These results can also be compared to a recent investigation 
by the Reference In put Parameter Library (RJPL) community 
dGorielv et alj|2009l) . making use of the code TALYS and HFB- 
14 fission barriers plus nuclear level densities obtained from a 
combinatorial approach based on single particle spectra from the 
corresponding Hartree-Fock-B ogolyubov calculations. The fis- 
sion barriers in dGorielv et al.ll2009l) were employed following 
the fission path via a Wentzel-Kramers-Brillouin phase integral, 
which should be superior to a Hill-Wheeler inverted parabola 
treatment. When examining their Fig. 8, which is based on pure 
predictions, it can be recognized that our results with HFB-14 
fission barriers employed via a double-humped fission barrier 



approach are similar or even closer to the experiment than the 
results given therein. As we make use of the same barriers, we 
relate the difference mainly to the different level density predic- 
tions (here a back-shifted Fermi gas with parameters obtained 
from a global mass model, there a combinatorial approach to sin- 
gle particle spectra from a microscopic HFB mass model). We 
draw the same conclusions from the comparison with their Fig. 
9 and our calculations with experimental barriers. While micro- 
scopic investigations should be in principle more advanced, the 
back-shifted Fermi gas approach based on global mass models 
seems still more robust in its predictive power. While renormal- 
izations of fission paths and level densities can improve this pic- 
ture (see their Figs. 10-12), this is only possible if experimental 
data are available. As our investigations are meant for astrophys- 
ical applications far from stability, where no experimental infor- 
mation is available, we come to the conclusion that our approach 
is well suited for this endeavor. It does, however, depend on the 
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quality of mass and fission barrier predictions, and this will be 
the focus of the present work. 

In astrophysical plasmas, reactions also occur on thermally 
excited target states. The stellar cross section <x* can then be 
defined as a sum of the cross sections o~ x for those excited 
states x with excitation energy E x and spin J x , weighted by the 
Boltzmann excitation probability 



Z x (27,+ po-V 
I,(2J x +l)e-l 



(2) 



The <x A were calculated in the same way as shown in Eq. ([TJ 
for the ground state, i.e. for x = 0. Only the stellar cross sections 
can be used to compute the appropriate astrophysical reaction 
rates. The reaction rate for a specific reaction at a given stel- 
lar temperature T was then determined by folding the stellar re- 
action cross section cr*(E) with a Maxwell-Boltzmann distribu- 
tionof relative velocities between projectiles and targets dFowlerl 
1974): 



<cr*v> = (ay)* 



1/2 

'V-i ! (kT) 



oo 

l — J a-*(E)Eexp(-^JdE. 



(3) 



Figures [2] and [3] show a typical comparison of present neu- 
tron capture rate calculations of 23 8 U and 242 Pu to expe rimentally 
based rates f rom JENDL-3 . 3 dNakagawa et al.l [2005b and th e 
predictions of lAikawa etal] (120051) : iGorielv et all (120081, 120091) . 
The agreement between rate predictions and data for other nu- 
clei in this mass range are of the same order. Note that these pre- 
dictions along the valley of stability are based on experimental 
masses, which leads to an average agreement with experimental 
cross sections within a factor of 1.5. The expected variation will 
be larger far from stability where theoretical mass and barrier 
predictions have to be utilized. 

3. Neutron-induced fission rates for a variety of 
mass models 



Early r-process calculations (Thielema nn et al.ll 19831) , which in- 
cluded fission, made use of the mass predictions by iHilf et all 
d 19761) an d the fission barri e rs of a macroscopic-microscopic 
model by iHoward & Mollerl (Il980h . For many years differen t 
authors used the fission barriers from IHoward & Mollerl (11980). 
as they were the only complete set of barriers available. More 
recently, renewed interest (and increased computing power) 
spurred a number of new calculations of large sets of barrier pre- 
dictions within various models, resulting on " average" in higher 
values of fission barriers than predicted by Howard & Mollen 
(1980). For a consistent treatment of nucleosynthesis, fission 
rates should be calculated with the neutron separation ener- 
gies, reaction Q-values and fission barrier h eights derived from 
the same mass mode l (see the discussion in ICowan et al.lll99ll 
iRauscher et al]|1994|) . 

As explained in Sect. |2 different mass- and fission barrier 
predictions were utilized to test the sensitivity stemming 
from different underlying models. The models used are: 
mass es taken from the Finite Range Droplet Model (FRDM) 
by iMoller et al.l (Q995), the Exte nded Thomas-Fermi with 
Strutinsky Integral (ETFSI) model bvlAboussir et all d 19951). an d 
the Thomas-Fermi (TF) model of iMvers & Swiateckil (1 1 996): 
fissio n barriers are taken from the ETF SI dMamdouh et al.l 
1200 lb and TF (IMvers & Swiateckil 1 19991) models. It should 
be emphasized that the ETFSI masses employed here 



(http://www-astro.ulb.ac.be/Nucdata/Masses/etfsi2-plain) 

are based on the force SkSC18 (IGorielvl 
2000), while the ETFSI fission barriers 
(http://www-astro.ulb.ac.be/Nucdata/Fisbar/fisbarr) were 
obtained with the force SkSC4 (Mamdou h et al.lll998l) . Thus, 
it is difficult to perform fully consistent calculations, and the 
quality of barrier heights is not really known especially in 
the region far from experimentally known nuclei, where the 
r-process proceeds. For this reason the choice of two different 
sets of fission barriers, TF and the ETFSI, permits to test the 
sensitivity range. When utilized together with mass predictions 
from the same models (see however the remarks in the previous 
paragraph), a reasonably consistent treatment is possible. In 
addition, we also chose to explore a combination of TF fission 
barriers and FRDM nuclear mass predictions and to gauge 
the resulting effect. The latter are close to the TF predictions, 
conta ining also the same shell corrections (My ers & Swiateckil 
1996) and have already been used in some astrophysical rate 
calculations for nuclei with Z < 84. 

Figure |4] shows the quantity of primary importance, Bf -S„, 
for the calculation of neutron-induced fission cross sections, 
based on the two model sets FRDM+TF and ETFSI. Bf - S n in- 
dicates the regions of nuclei where neutron induced fission rates 
can be high enough (Bf - S„ < 0) and are important for the r- 
process nucleosynthesis (see also the discussion of Fig. 6). We 
see that both sets display a quite different behavior and note that 
Bf-S„ is generally larger for ETFSI than FRDM+TF. 

Our fission cross section calculations were performed 
within the framework of a double-hump fission barrier (per- 
mitting the existence of a double or single hump barrier). The 
he ights of both barr i ers f were predi cted in the calculation s 
of IHoward & Mollerl (Il980h (HM) and lMamdouh et al.l d2001h . 
However, the TF model predicts only one of the fission barriers. 
In order to employ this model in our calculations we assumed 
that the predicted fission barrier is the higher of bo th barriers 
and de rived the lower one as described below and in lPanov et al.l 
(12001 . In order to derive the height of the lower barrier, we 
compared two methods: (i) keeping the difference in height of 
the first and second barrier of the original HM barriers and (ii) 
keeping the same height ratio, i.e. relative height, of the barri- 
ers. The resulting cross sections were not very different because 
the heights of the lower barriers calculated in these two ways 
differ only by a few percent for the majority of cases. Only in 
a small number of cases the cross sections differ by more than 
10% (but the largest deviations do not exceed a factor of two). 
For the rates shown here we chose to use the difference in height 
of the first and second barrier of HM to predict the lower barrier 
for the TF approach. In this manner, the new fission rate calcula- 
tions were extended to th e region above charge Z=100, whereas 
IHoward & Mollerl (Il980l) presented results only for Z < 100. 
This allows us to perform r-process calculations in extremely 
neutron-rich environments as well as to study superheavy ele- 
ment production in rapid neutron capture nucleosynthesis pro- 
cesses. 

An extended comparison to evaluated neut ron-induced fis- 
sion cross sections, based on evaluated data (Nakagawa et al.l 
2005) for the trans-lead region is shown in Fig. [5] The left panel 
displays the ratio of predicted to evaluated cross section when 
using experimentally known fission barriers (circles). The agree- 
ment of the calculated rates based on the experimental values of 
fission barriers is quite good with the majority of ratios, and is 
within factors of two to three. Some outliers of up to a factor of 
ten are observed in a few cases, but the experimental accuracy 
of the barrier determination is not known for these cases. The 
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Fig. 2. Comparison of (n, y)-rates (integrated over Maxwell-Boltzmann distributions of targets and projectiles for the displayed tem- 
peratures) from present calculations and other existing pre dictions for the target nuclei 238 U (left) and 242 Pu (right) with experiment. 
The symbols are chose n as follows: lAikawa et alj d2005l) (blue dash-do t line). iGorielv et all {2008) (thin blue line), JENDL-3.3 
dNakagawa et al.ll2005l) (crosses), new Talys-based (iGorielv et alJ 120091) predictions (green dash-dot line). Our calculations are 
shown for reaction rates including only the g round state (GS) or a thermally populated target (GS+tps). Only the GS rates can be 
compared to experimental data. The results of lAikawa et al.1 (120051) : IGorielv et all d2008l) correspond to GS + tps conditions. 




Fig. 3. Comparison of (n,/)-rates (integrated over Maxwell-Boltzmann distributions of targets and projectiles for the displayed 
temperatures) from present calculations and ot her existing p r edicti ons for the target nuclei 238 U (left) and 242 Pu (right) with the 
experiment. The symbols are chosen as follows: IGorielv et al.l d2009) (green two-dot-dash line), JENDL-3.3 (Nakagawa e t al.ll2005l) 
(crosses). Our calculations are shown for reaction rates including only the ground state (GS ) or a thermally populated target (GS+tps) 
for different nuclear data predictions (HFB or ETFSI). The results of IGorielv et al.l (120091) correspond to GS + tps conditions. 



middle panel shows the difference between the calculated ones 
when the HFB mass and fission barrier predictions were used: in 
this work (triangles down) and by Talys (crosses). Existing dif- 
ference in two Hauser-Feshbach calculations can have emerged 
from the differences in fission barriers values (due to renormal- 
ization by IGorielv e t al. (2009)), different level density used etc. 
In the right panel of Fig. [5] the calculated cross sections made 
use of mass and fission barri er predictions from TF and ETFSI 
as well as new predictions by Mfall er et"a"D (120091) . 

Contrary to the comparison with the left panel (circles), 
when using only experimental barriers we found that the differ- 
ence between calculated cross sections and measured ones can 
be many orders of magnitude (up to a factor 10 4 ). It is interest- 



ing to note that a weak systematic can be seen. Cross sections 
calculated with ETFSI as well as HFB fission barrier predictions 
show a large scatter above and below the measured values for 
the lighter end of the plotted mass range, whereas they tend to 
underpredict the fission cross sections f or heavy ma s ses. Th e cal- 
culations using the recent barriers by Moll er et"a"fl d2009) seem 
to show a similar behavior with a smaller scatter, but they are 
currently available only for a more limited number of nuclei. 
The results obtained with the TF barriers exhibit a different pat- 
tern, the predicted cross sections agree well for the lighter mass 
nuclei, but seem to be systematically too large for the heavier 
ones. 
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Fig. 4. Contour plot of the quantity Bf - S„ (in MeV) for two choices of models for masses and fission barriers: FRDM masses 
plus TF fission barriers (left) and ETFSI-predictions (right). Also indicated are known superheavy nuclei (boxes) and the drip-line 
(line), where Bf—S n determines whether the r-process continues towards heavier nuclei or fission cycling to lighter fission products 
occurs. 
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Fig. 5. Calculated maxwellian averaged neutron-induced fission cross sections in comparison to evaluated experimental fission 
cross sections at 30 keV. Left panel: experimental fission barriers were utilized in the calculations (circles) of m axwellian aver- 
aged cross-sections in comparison to experimental values for 58 isotopes from Th to Fm taken from JENDL-3.3 (Nakagawa et al. 
2005); central panel: the same ratio of calculated and evaluated experi mental values, bu t emplo ying theoretical fission barrier pre- 
dictions from different sou rces: the Extended Thomas-Fermi m odel by [M amdou het al.l d200ll) (diamonds), the Thomas-Fermi by 
Mvers & Swiatecki ( 1999) (squares), and recent predictions from Moller et al.l (l2009l) (green triangles up). When not experimentally 
known, consistent nuclear masses were utilized from the corresponding mod el predictions. Right panel: comparison of maxwellian 
averaged cross-sections from our calculations (triangles) and from the Talys (Gori elv et al.l l2009) calculations (crosses), both utiliz- 
ing HFB predictions (with BSkl4 Skyrme force) of masses and fission barriers. 

Table 1. Neutron separation energy S„ and fission barrier predictions for nucleus 262 U formed after neutron capture by 261 U. 



Mass and Bf models 


ETFSI 


TF 


TF FRDM 


TF ETFSI 


HM 




Bfl Sf2 S „ 


Bfl Bf2 S n 


Bf\ Bfi S „ 


Bfl Bf2 S „ 


Bfl Bf2 S n 


calculated values, in MeV 


3.9 5.30 4.46 


1.20 4.56 4.05 


1.20 4.56 3.81 


1.20 4.56 4.46 


3.36 4.14 



The above comparison underlines the considerable varia- 
tions still inherent in fission barrier predictions. However, we 
sugge st that by comparing TF with ETFSI (and with Moller et~a"fl 
2009) predictions, the relevant variation range can be estimated. 
Theoretical cross sections depend strongly on the fission barri- 
ers, and a high accuracy for their values is required. Because of 
the impact of the fission barrier uncertainties, nucleosynthesis 
studies at present should explore a variety of barrier sets, while 
waiting for further measurements and improved predictions. For 



this reason we compute and compare below rates for different 
sets of nuclear properties and also provide tables and fits of the 
rates for all these cases. 

Figure [6] shows fission cross sections (left panel) and rates 
(right panel) calculated by combining different sets for predic- 
tions of masses and fission barriers. The arrows in the left plot 
show the difference between the fission barrier height and the 
neutron separation energy Bf — S„ given by predictions of the 
TF model (red arrow at the top of the left panel) and the ETFSI 
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model (dashed arrow at the bottom). The exact values can be 
found in Table Q] The cross sections (and thus the rates) depend 
essentially on the available energy Bj - S„ in the fission chan- 
nel, minor dependencies on fission barriers heights and neutron 
separation energies, individually, are due to the competition with 
the (n, y)-channel. 

It can be seen from Fig. [6] that the neutron-induced fission 
cross sections (as well as the rates) increase with decreasing fis- 
sion barriers and that cr™ > crJJ > cr^J FS1 . The difference at 
low energies is due to the different mass predictions used (com- 
paring calculations with TF fission barriers but different neutron 
separation energies S n ). The small decrease in S n , when predic- 
tions for masses and fission barriers (based on the two sets of 
input data) are changed from TF+TF to FRDM+TF, results in a 
decrease of the neutron-induced fission cross sections (accom- 
panied by decreasing (n, y)-cross sections). The same influence 
is illustrated by the cases where sets of consistent determinations 
for S„ and Bj are replaced by sets from different mass predic- 
tions. For example, different predictions of neutron separation 
energies were used for the same TF - fission barrier predictions. 
Therefore, cross sections for low energies differ significantly, 
but for higher E the difference in cross section becomes much 
smaller for the same fission barriers (here we used TF barriers). 
The temperature averaged rates (Fig. [6] right panel) show the 
same dependence. 

As expected, the rate sets calculated on the basis of FRDM 
masses and the TF model for masses and fission barriers are 
quite comparable (Fig. |7). On the other hand, when the sets in- 
clude ETFSI vs. FRDM+TF or TF comparisons, the results can 
differ by up to eight orders of magnitude. The most extreme dif- 
ference of ETFSI-based and other rates was obtained for nuclei 
with neutron numbers close to 184, for which the ETFSI model 
predicts very high fission barriers, leading to small fission rates. 
The difference between rates, calculated on the basis of ETFSI 
and HFB mass and fission barriers predictions is less than be- 
tween ETFSI and TF, especially for higher Tg. For smaller Tg 
the we can see that difference of rate values has the opposite sign 
for regions with A w 240 - 280 and A > 280. This agrees with 
Fig- El Some of the combinations employed in Fig.|7]are shown 
to underline the huge problems which can arise when not using 
consistent data. The extrapolation of rate calculations to regions 
of very exotic nuclei is a hard task, and only further investiga- 
tions can answer which kind of prediction is more preferable. At 
the moment the only choice is to test all available predictions in 
r-process calculations that are compared to astronomical abun- 
dance observations. 



4. Computed rate sets and mass ranges 

As we discussed in the introduction, our aim was to prepare 
sets of neutron-induced reaction rates, suited for calculations 
of the r-process up to very high atomic masses. In this sense, 
our rates extend the previousl y published ones for Z < 84 
dRauscher & Thiele mann 200l|) up to the region 84 < Z < Z max 
(Z max depends on available nuclear input and varies from 102 for 
HFB-predictions to 1 18 for FRDM data). The exact range of nu- 
clei for the different choices of mass and fission barrier input is 
given in Table [2] 

Our calculations include all outgoing channels and give 
simultaneously predictions for neutron-induced fission, (n,y)-, 
(n,p)-, and (n,a) rates. However, here we only provide tables 
and fits for the neutron-induced fission and neutron capture rates. 
Due to their small values for neutron-rich nuclei, the rates of the 
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other calculated reaction types are not important in astrophysical 
applications. 

The format of the tables is explained in Appendix A. The 
full electronic versions of the tables available on-line at the CDS 
include all rates for all mass predictions (Tables 3-6). The iso- 
tope and element ranges for which rates as well as rate fits are 
available are given in Table [2] for the FRDM, TF, ETFSI and 
HFB mass predictions. This amounts to 2151 (ETFSI), 2637 
(TF), 2400 (FRDM-masses, TF-barriers) and 1323 (HFB) in- 
volved nuclei. The partition functions for all isotopes are given 
on a grid of 24 temperatures: T 9 = 0.1, 0.15, 0.2, 0.3, 0.4, 0.5, 
0.6, 0.7, 0.8, 0.9, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 6.0, 7.0, 

8.0, 9.0, 10.0 and can be found in CDS's Tables 3-6 as well. We 
also provide the fit coefficient needed to compute the photodis- 
integration rate (see Sect. 14. II ). 

4.1. The fits for neutron-induced and reverse rates 

Reaction rates have been calculated on a grid of 24 tempera- 
tures: r 9 =0.1, 0.15, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.5, 
2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0, same as 
for partition functions. These rates include the thermal modi- 
fication in the stellar plasma, i.e. consider reactions from ex- 
cited states in the target. For easy application in astrophysi- 
cal investigations, these stellar rates were fitted with the same 
REA CLIB parameterization a s used for other reaction types ear- 
lier dRauscher & Thielem ann 200Cj): 

Na (o~v)* — exp (ao + a\T^ x + a^T^^ + a^T^ 3 + a\Tg 

+a 5 T^ 3 +a 6 lnr 9 ) , (4) 

with the seven open parameters ao - a& and the stellar tempera- 
ture Tg given in 10 9 K. This parameterization proved to be flex- 
ible enough to accommodate the different temperature depen- 
dencies of the various reaction types across the fitted tempera- 
ture range of 0.01 < Tg < 10. Paramet e rizatio ns of the present 
rates in the form used in Holme s et al ] (Il976h can be found in 
Appendices. The best fit was obtained b y minimizing the d evia- 
tion £ (Eq. Q) using the FUMILI code iDvmov et alhOOOl) . 

The flexibility of the fitting function makes it prone to nu- 
merical problems outside the calculated range at low tempera- 
tures. In some cases they tend to diverge strongly. This difficulty 
can be avoided by additionally providing fit data at low tem- 
peratures to the calculated values by appropriately extrapolating 
the rates to lower temperatures. However, it has to be empha- 
sized that the considered parameterization is only valid within 
the temperature range of 0.01 < Tg < 10.0, although many fits 
will show a "proper" behavior down to lower temperature. 

As a measure of the accuracy of a given fit, the quantity £ 
(marked in tables and figures as Dev) is used. It is defined by 




with r being the original rate value as calculated at each of the 24 
temperatures Tg- 0.1, 0.15 ... 10.0, and fi is the rate calculated 
from the fit at these temperatures. A small value of £ indicates 
an accurate fit over the entire temperature range. Higher values 
of ( are mainly caused by deviations at low temperatures, where 
rates are slow and a larger deviation is permissible. For the ma- 
jority of nuclei the value of £ is less than 1 and lies in the range 
0.1 - 10 (see Fig.[9j». We should mention here that t he modified 
approximation formula by Cvbur t & Davids! d2008l) can proba- 
bly give the better fit. Its accuracy can be high and should be 
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Fig. 6. Dependence of neutron-induced fission cross sections cr n f(E) (left) and rates A n f = Na (<xv)(right) on mass- and fission 



barrier predictions for U. The sources for fission barriers Bf and neutron separation energies S„ (Howard and Mtiller, Thomas- 
Fermi, extended Thomas-Fermi, finite range droplet model) are indicated in the panel. Arrows show the difference between fission 
barrier and neutron separation energy Bf — S„ for ETFSI (dashed line) and TF (full line) predictions. 
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Fig. 7. Ratio R = of reaction rates for different fission barrier and mass predictions at temperatures Tq = 0.3 (left) and 

Tg = 10 (right). The indexes I,J define the mass/barrier predictions utilized in the calculations: ETFSI (I,J=1), Thomas-Fermi (I,J=2) 
or Thomas-Fermi for fission-barriers and FRDM for mass-predictions (I,J=3), Hartree-Fock-Bogolyubov (HFB-14) predictions 
(I,J=4). As the Thomas-Fermi model and the Finite Range Droplet Model give similar mass predictions, the deviation for the 
combination 2/3 is the smallest (and also decreases with temperature). However, the use of TF (or FRDM) vs ETFSI (or HFB-14) 
predictions changes the rates drastically and thus the full range of possible magnitudes shown in Fig.|4]is obtained. 



applied for reactions, whose rates can be calculated with rather 
high accuracy. For our case fitting the average accuracy is not 
bad (see Fig. [8j, and is much better than accuracy of calcula- 
tions of neutron-induced rates for very neutron rich nuclei. We 
applied the approximation Eq. (0J, used earlier in a number of 
previous predictions of neutron rates. 

The temperature dependence of the rate can be one of two 
types, as illustrated in Fig. [8] These types of behavior can be 
understood when one recalls the discussion of Fig. [6] A fission 
transmission coefficient which is constant or slowly varying as a 
function of energy leads to an (n, /)-cross section which (similar 
to a pure neutron capture) shows a 1 / Ve dependence, if s-wave 
dominated. Averaging such a cross section over a Maxwell- 
Boltzmann distribution yields a constant rate. This situation oc- 
curs for example when the neutron bombarding energy leads to a 
compound nucleus energy above the lower and below the higher 
barrier of a double-hump fission barrier. Then the lower barrier 



is open and the penetration through the remaining higher barrier 
is close to constant, but the height of the higher barrier deter- 
mines the size of the (n,f) cross section. This behavior is seen 
in Fig. [6] below about 0.1 MeV (left panel). In the right panel 
the corresponding rate is shown and seen to be close to constant 
below about Tc, — 1 . 

In contrast, the opening of a fission barrier as a function of 
energy, i.e. an exponentially increasing transmission coefficient 
close to the barrier energy, leads to a sudden rise of the cross 
section, and consequently also of the rate as a function of tem- 
perature. We see this behavior for bombarding energies on the 
order of Bf - S „ in the left panel of Fig. |6]and a similar behavior 
for the rates in the right panel. The size of this change from an 
almost constant to a steeply rising rate is a function of the barrier 
height. The right panel of Fig.[8]is a very representative exam- 
ple of such a case (almost constant rate at low temperatures and 
a steep rise by orders of magnitude at a critical temperature). 
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Table 2. Ranges of isotopes for which we calculated rates based on input from the models TF, FRDM, ETFSI and HFB. Given are 
the charge number Z and the lower and upper limits A m ; n and A max of the neutron number of the targets in the isotopic chain. 





TF-shc 


FRDM 


ETFSI 


HFB 




TF-shc 


FRDM 


ETFSI 


z 


4 . 


4 

'Mnax 


4 . 

"min 


4 


4 . 

^min 


4 

^max 


4 . 


4 

"max 


Z 


4 . 

^Mnin 


4 


4 . 


4 

^max 


4 . 


4 

^max 


84* 


180 


269 


180 


269 


180 


267 


193 


269 


102 


233 


331 


233 


331 


233 


331 


85* 


185 


270 


185 


270 


185 


270 


193 


270 


103 


240 


335 


240 


334 


240 


328 


86* 


187 


269 


187 


269 


187 


269 


196 


269 


104 


239 


337 


239 


337 


239 


324 


87* 


190 


280 


190 


280 


190 


280 


200 


278 


105 


242 


337 


242 


331 


242 


330 


88* 


193 


283 


193 


283 


193 


279 


202 


279 


106 


245 


337 


245 


331 


247 


337 


89* 


196 


288 


196 


288 


196 


284 


206 


285 


107 


248 


337 


248 


332 


243 


337 


90 


199 


293 


199 


293 


199 


287 


220 


287 


108 


251 


337 


251 


327 


245 


337 


91 


212 


296 


212 


296 


212 


288 


229 


289 


109 


254 


337 


254 


327 


247 


337 


92 


204 


299 


204 


299 


204 


291 


230 


291 


110 


257 


337 


257 


339 


249 


337 


93 


220 


302 


220 


302 


220 


293 


232 


293 


111 


260 


337 


260 


331 


251 


329 


94 


210 


305 


210 


305 


210 


295 


234 


295 


112 


263 


337 


263 


332 


253 


329 


95 


215 


309 


215 


308 


215 


297 


238 


297 


113 


267 


337 


267 


327 


255 


329 


96 


216 


312 


216 


311 


216 


300 


240 


300 


114 


270 


337 


270 


327 


257 


329 


97 


225 


315 


225 


314 


225 


303 


243 


315 


115 


273 


337 


273 


337 


275 


329 


98 


222 


319 


222 


319 


222 


317 


245 


319 


11 6 


276 


337 


276 


337 






99 


225 


322 


225 


322 


225 


320 


250 


322 


117 


279 


337 


279 


337 






100 


227 


325 


227 


325 


227 


325 


251 


325 


118 


282 


337 


283 


337 






101 


239 


328 


239 


328 


239 


326 


254 


328 


119 















* for these chemical elements there are no HFB-predictions of fission barriers 




log 10 (T 9 ) log 10 (T 9 ) 

Fig. 8. Representative fits of neutron-induced fission rates ( 273 Cf and 291 Cm) for two types of temperature dependences. The general 
aim is to attain sufficient accuracy with one set of seven fitting coefficients (see Eq. (4)), as shown in the left curve. In some cases 
the superposition of two such sets is necessary (right panel) to correctly reproduce the low- as well as high-temperature behavior. 
A correct extrapolation to low temperatures (0.01 < Tg < 0.1) is important to avoid unphysical abundance changes (see also the 
explanation of tables). 



On the other hand, when the compound energy is close to the 
barrier height already for small bombarding energies, the cross 
section and rate are already large at small energies (tempera- 
tures). Such an example can be seen in the left panel of Fig. [8] 
reflecting a double-hump behavior with small energy differences 
between the lower and higher barrier. For a discontinuity in the 
7Vdependence (right panel), the fit was performed as a sum of 
two contributions and is given by two lines in the Table IA.2I of 
the paper and CDS's Tables 7-18. 

For all cases it is recommended to use the fits only down to 
the temperature Tg > 0.01. Moreover, close to the drip-line, the 
statistical model may not be applicable for reactions with low 
Q value, even above that temperature. Although the fit may be 
good, the user should be aware of that possible complication. 



4.2. Photodisintegration rates and partition functions 

For a full implementation of the neutron captures in a reaction 
network, the inverse photodisintegration rates also have to be 
known. The photodisintegration rates (and their fits) are not ex- 
plicitly given in the tables, but can be computed from the in- 
formation contained therein. To calculate the reverse rate of the 
reaction B(n,y)D, i.e. the reaction D(y,n)B, the seven parameters 
a T Q V - flg ev are determined as follows: 



= qq + In 



9.8685 x 10' 



,\ D - ! \" 27 B + 1 

2/n + 1 



(2/„ + 1) 



= a x - 11.60455,, 
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Fig. 9. The accuracy of fits is defined by a mean-square error (see Eq. (|5]l), shown here for (n, y)-rate fits (left panel) and neutron- 
induced fission rate fits (right panel). We see that the upper boundary of errors is in the range of 10%. This should be compared to 
the variations of these calculations (see Figs. Q] and |5J. 



a™ = a 6 + 1.5 , 

with Ad the mass number of nucleus D, / n the spin of neutron 
and Jb, Jd the ground state spins of nuclei B and D, respectively. 
These parameters for the reverse reaction are also given in the 
tables (see explanations in the appendices). 

It is important to note that the value computed by applying 
Eq. (H]i with the above coefficients has to be multiplied by the 
ratio of the partition functions for the residual and target nucleus 
Gb/Gd to obtain the actual photodisintegration rate. Examples 
are shown in Appendix [ A| A s it was shown in detail earlier 
dRauscher & Th ielemann 2000), the temperature-dependent par- 
tition function G(T*) n ormalized to the grou nd state spin J° of a 
nucleus is defined as in lFowler et al.l d 19671) 

(2J° + V>G(T*) = Tl^LoW + l)e- E " /kT " (7) 

+ / E7^(2^ + l)e- e/fcr V(e,/^)d<: , 

with p being the level density and fi,„ the last included exper- 
imentally known state. For the temperature range considered 
here, the maximum energy £™ ax above which there are no more 
significant contribut ions to the partition function is on the order 
of 20 - 30 MeV dRauscher & Thielemannll2000l) . The tempera- 
ture dependent partition functions are available at the CDS, as 
well as the fit coefficients in Eq. (0]i (Tables 7-18) for the (n, y), 
(y, n) and (n,f) rates. 

This subsection discussed photo-induced reaction rates as 
inverse reactions of neutron capture rates, which are presented 
in full detail in the present publication. Photo-induced reactions 
can also lead to fission when the Planck distribution of photons 
provides a significant fraction of photons with energies above the 
fission barrier. This reaction channel is not discussed here but its 
possible influence in r-process environments will be analyzed in 
a future investigation. 



5. Fission fragment distributions 

A proper inclusion of fission in r-process calculations also re- 
quires the knowledge of the resulting distribution of fission frag- 
ments, which have to be entered as reaction products. This by it- 
self requires a major effort and will be presented with a thorough 
description of the treatment plus detailed fission yield distribu- 
tions in a forthcoming paper. However, at the end of the present 
investigation we want to give a short outlook on how this topic 
will be approached. 

In principle, such distributions are dependent on the exci- 
tation energy of the compound nucleus and thus would lead 
to a changing yield distribution for each bombarding energy. 
However, we found that the distributions vary smoothly and 
slowly as a function of the excitation energy. Thus, for the neu- 
tron energy range in astrophysical applications, the yield dis- 
tribution at the neutron separation energy (i.e., for a vanishing 
neutron bombarding energy) is a very good approximation. This 
would permit us to multiply the fitted (n,f)-rates from the pre- 
vious section with a static distribution of yields for all tempera- 
tures. 

Here we only briefly show a few examples of such fis- 
sion distributions. Properties of fission fragments, i.e., masses, 
atomic numbers, excitation and kinetic energies, were calculated 
based on the macro-microscopic approach (similar to the FRDM 
model) and the separability o f compound-nucleus and fragment 
prope rties on the fission path ( Sch midt et al.' 2008; Wilkins et al.l 
1976). The original technical de scription of the frag ment- 
formation model wa s published in iBenlliure et al.l d!998l) and 
iRruglov et al.l d2002l) . In the calculations shown here we used 
an updated description that will be the subject of a forthcoming 
publication. 

In the model it is assumed that the different ways of split- 
ting up the total mass are basically determined by the number 
of available transition states above the potential energy surface 
behind the outer saddle point. The macroscopic properties of the 
potential-energy landscape of the fissioning system are attributed 
to the strongly deformed fissioning system, which are deduced 
from mass distributions at high excitation energies dRusanovl 
fl997h and Langevin calculations (Nadto chv et alj2005l) . The mi- 
croscopic properties of the potential-energy landscape of the fis- 
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sioning system are given by the qualitative features of the shell 
structure in the nascent fragments. They ar e determined from 
the observed features of the fission chan nels dBrosa et al.ll 1990) 
according to the procedure described by Schmidt et al. (2008). 
The dynamics of the fission process responsible for the frag- 
ment formation was considered in an approximate way: It was 
assumed that the phase space near the outer saddle point de- 
termines the mass asymmetry of the system, which is more or 
less frozen during the descent to scission. On the other hand, it 
was also assumed that the N/Z collective degree of freedom is 
determined near the scission point. The excitation energies of 
the created fragments were calculated from the available excita- 
tion energy at the scission point and the deformation energies of 
the fragments at scission. The deformation energies of the frag- 
ments were assumed to be specific to the individual fission chan- 
nels. They were deduc ed from experimental data dWahll 1988: 
lBockstiegelet"ai] |2008. and references therein) on total kinetic 
energies and neutron yields. Kinetic energies were then calcu- 
lated applying the energy conservation law. Finally, the two ex- 
cited fission fragments are subject to particle (mostly neutrons) 
and y ray emission until they reach their ground state configura- 
tions. The de-excitation process was de scribed in the fram ework 
of the statistical model as described in iKelic et al.1 d2008l) . 

The left panel of Fig. [10]compares the mass distributions re- 
sultin g from such an a pproach with the experimental data for 
238 U dNagy et aLlll978l) . We also show the yields c omputed us- 
ing th e e mpirical parameterizations de veloped by iPanov et"aTI 
( 2008f) and lKodama & Takahashild 19751) . which previously were 
used in r-process calculations. In order to give an impression of 
the impact of the results for r-process nuclei, the two right panels 
of Fig. [10] show the mass distributions resulting from the fission 
of 282 Cm and 292 Cf. For these nuclei, clear differences appear 
between the phenomenolog ical approaches dPanov et al.l [2008: 
iKodama & Takahashil fl975) and the macroscopic-microscopic 
approach described above. For application in r-process sim- 
ulations the (n,/)-reaction rates (which are the focus of the 
present paper) will have to be multiplied by the displayed per- 
centages for o btaining a produ c tion rate of a specific frag- 
ment nucleus (IKelic et al.l l2005t: iMartfnez-Pinedo et al] l2007t 
IPanov et al.ll2008l) . 

6. Summary and conclusions 

We provide predictions of neutron-induced fission rates and 
(n,y)-rates for a wide range of astrophysical temperatures 
(10 8 < T(K) < 10 10 ) and targets (proton- to neutron- drip- 
line for 84 < Z < 118, i.e. from Po to Uuo ) in the 
framework of the Wolfenstein-Hauser-Feshbach model, mak- 
of a variety of different mass and fission barrier 



use 



predictions (|M\ 


/ers & Swiateckil 1 19991 Mamdouh et al.| 20011 


Aboussir et al. 


19951 iMoller et al. 11995b iHoward & Moller 


1980: iGorielv et al. |2009|). The astrophysical (stellar) reaction 



iRauscher & Thielem ann 2000) in the common REACLIB seven 
parameter form, and these parameters are also tabulated. This 
provides the basis for r-process nucleosynthesis calculations 
where the abundance predictions for the highest mass numbers 
as well as the effect of fission cycling are strongly dependent on 
the interplay of neutron capture and fission. 

In order to give an impression of the reliability of the results, 
we compared them with experiment and with available indepen- 
dent predictions before exploring the currently unreachable re- 
gions of the nuclear chart with a variety of theoretical predic- 
tions for nuclear masses and fission barriers (FRDM, ETFSI, TF, 



HFB). An extended comparison of neutron-induced fission rates 
with experiment and with available independent predictions was 
done. The dependence of rates on nuclear input data, most of all 
fission barriers, is high. Astrophysical nucleosynthesis yield pre- 
dictions, especially in the transuranium region, should take into 
account these large differences in order to explore the variations 
involved. For this reason extended tables for neutron-induced fis- 
sion rates as well neutron capture rates are presented for different 
mass and fission barrier predictions in fitted form for nucleosyn- 
thesis calculations. Their structure is given in the Appendix A 
(note that the full rate and fit tables are available at the CDS). 
Given that fission predictions far from stability have not been 
tested yet, and even close to stability none of the existing models 
has yet been proven to be superior (see Fig. |3j. Nucleosynthesis 
calculations should probably continue to use a variety of these 
models. A further requirement for nucleosynthesis modeling in 
the region of fissioning nuclei is the knowledge of the mass dis- 
tribution of fission products. This work is in progress (see Sect. 
0. 
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Appendix A: Explanation of the tables and 
examples of how to use them 

This section is intended to help with interpreting the informa- 
tion given in the electronic tables. We also give examples for 
calculating the reaction rate for a given reaction and its inverse 
reaction. 

Table IA.1I shows which information is contained in the rate 
Tables 3-6 for different mass predictions, available at the CDS 
(Table |A~T1 illustrated some lines from CDS Table 3: rates on 
the basis of ETFSI-predictions). The organization of Tables 3-6 
with calculated rates for 24 different values of Tg is extremely 
simple - all data are in eight columns. The columns give, in this 
order, the target element, the atomic mass number A of the tar- 
get, the temperature in Tg, the partition function of the target, the 
neutron-induced fission rate for the ground state, the neutron- 
induced fission rate with thermally populated target states, the 
neutron capture rate for the ground state, and the neutron cap- 
ture rate with thermally populated target states (units of all rates: 
cm" 3 mofe _1 c _1 ). 

The rates in Tables 3-6 were calculated on the basis 
of different mass and fission barrier predictions: ETFSI, TF, 
FRDM(masses)+TF(barriers) and HFB respectively (see main 
text of the paper for details). 

The entries are denoted as follows: 

mother mother nucleus (target), 
Tij stellar temperature in 10 9 K, 

p.f. nuclear partition function, 

(n,g) (n,y)-rate, N A (cr„ y v) 

(n,g)* (n,y)-rate, Na (cr„ r v^, with thermally populated tar- 
get levels, 
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Fig. 10. The final mass distributions of fission fragments for the compo und nuclei (after neutron capture) 238 U, 282 Cm and 292 Cf. 
The distributions were computed with the ABLA code (Kelic et al. 2008) as descr ibed in the text. In addition, we show the yields 
computed with the phenomenological parameterizations of Panov et al. (2008]) and lKodama & Takahashil (1 19751) . as well as experi- 
mental data (crosses). 



Table A.l. Example of format of rates presented in the Tables 3-6 , available at the CDS: reaction rates Na (o~v) and partition 
functions show a subset of calculation based on ETFSI predictions. 



mother(A,Z) 


T 9 


p.f. 


(n,f) 


(n,f)* 


(n,g) 


(n,g)* 


Cf 


273 


0.10 


1.00D+00 


5.07D+08 


5.90D+08 


2.89D+04 


3.32D+04 


Cf 


273 


0.15 


1.01D+00 


6.59D+08 


8.13D+08 


3.71D+04 


4.50D+04 


Cf 


273 


0.20 


1.02D+00 


7.81D+08 


1.01D+09 


4.35D+04 


5.50D+04 


Cf 


273 


0.30 


1.08D+00 


9.70D+08 


1.33D+09 


5.28D+04 


6.99D+04 


Cf 


273 


7.00 


7.86D+04 


1.89D+09 


1.57D+09 


5.70D+04 


9.86D+01 


Cf 


273 


8.00 


9.13D+05 


1.93D+09 


1.53D+09 


5.41D+04 


1.47D+01 


Cf 


273 


9.00 


1.04D+07 


1.98D+09 


1.60D+09 


5.14D+04 


2.20D+00 


Cf 


273 


10.00 


1.03D+08 


2.03D+09 


1.97D+09 


4.89D+04 


3.69D-01 


Cf 


274 


7.00 


3.87E+05 


1.75E+08 


1.68E+09 


4.04E+06 


5.85E+02 



(n,f) (n,fission)-rate, Na \c n fvj 

(n,f)* (n,fission)-rate, Na \ cr nfv\ with thermally popu- 
lated target levels. 

Table IA.2I shows which information is contained in the fit 
Tables 7-18, available at the CDS. The following information is 
provided: 

mother reaction target 
Dev fit accuracy £ (Eq. Q) 

flo. . . fl6 seven fit parameters for the forward rate 
a™. . . aj. ev seven fit parameters for the reverse rate fit (see 
Sect.O 

The fitting coefficients for the (neutron,y)-, (y,n)-, and 
neutron-induced fission rates with different mass and fission- 
barrier predictions are placed in the Tables 7-18. The columns 
in Tables 7-18 are organized as follows: target element, atomic 
mass number A of the target, target charge number Z, the number 
of fitting curves seven coefficients of the forward reaction a,-, 
and the mean square error. 

A value i/n-0 means that there is only one seven-parameter 
set to fit the rate. Values //,,> 1 give the number of parameter sets 
which have to be added up to yield the final rate, i.e. the rate r 



is calculated as r — 2, ri, with each r, computed from the z'-th 
parameter set and using Eq. (0]). 

The examples of Tables 3-6 are given for 1 isotope for 24 
values of Tg, and for Tables 7-18 - for 10 isotopes. Note that the 
tables of CDS show the values with an accuracy of seven digits. 

Below we give two examples for calculating a rate at Tg = 
7.0 with the fit parameters listed in the tables. 

The first example is the reaction 273 Cf(n,/). In Table A. 2 one 
finds the parameters ao = 38.72, a\ = -0.057, = 8.78, 03 

27.85, a A = 1.90, a 5 = -0.097, a 6 = 10.45. With the help of 
Eq. one calculates Na = 1.574 X 10 9 cm 3 s~ 1 mole~ 1 or 

U>g m (N A (ov)} u ) = 9.20 at T 9 = 7.0. 

The second example is for the rates of the capture 
reaction 273 Cf(n,y) 274 Cf and its reverse reaction. Similar to the 
above example, using the parameters from Table A. 2 and Eq. 
(0]l the capture rate is easily found to be N A (o~v)j it = 120 
cm 3 s _1 mole _1 . With the reverse parameters the first value in the 
determination of the reverse rate is found to be A 7 ' =5. Ox 10 9 s _1 
at 79 = 7.0. In order to obtain the actual value of the reverse 
rate, one first has to determine the ratio of the partition functions 
G273 Cf /G274 Cf = 4.58 x 10 5 /3.87 x 10 5 = 1.18 (see Table A.l and 
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Table A.2. Parameterization of the (n, y)-, (y, n)- and (n,f)-rates, available at the CDS Tables 6-18. 



Example of Table 7 — (n, y)-rate fits on the basis of ETFSI mass-model predictions. 



mother(A,Z) 


'fit 


ao 


a\ 


a 2 


a 3 


a 4 


a 5 


a 6 


Dev 


cf 268 98 





-1.45E+01 


0.00E+00 


-2.94E+00 


3.00E+01 


-4.69E+00 


1.40E-01 


-6.04E+00 


3.9E-02 


cf 269 98 





2.05E+01 


1.43E-02 


-6.95E-01 


-3.73E+00 


1.60E-01 


-3.13E-01 


7.64E-01 


5.0E-02 


cf 270 98 





-2.03E+01 


0.00E+00 


-3.86E+00 


3.56E+01 


-5.02E+00 


1.90E-01 


-7.77E+00 


1.9E-02 


cf 271 98 





-6.15E+00 


1.68E-02 


-3.20E+00 


1.94E+01 


-2.75E+00 


-3.24E-03 


-4.49E+00 


1.1E-02 


cf 272 98 





-1.69E+01 


n OOF +00 


-4.68E+00 


4.14E+01 


-5.82E+00 


1.64E-01 




3.6E-02 


cf 273 98 





2.62E+00 


1.71E-02 


-2.67E+00 


1.32E+01 


-1.60E+00 


-1.64E-01 


-3.27E+00 


2.2E-02 


cf 274 98 





-1.87E+01 


00E+00 


-5.35E+00 


4.34E+01 


-4.94E+00 


2.04E-02 


-1.03E+01 


4.9E-02 


cf 275 98 





2.43E+00 


0.00E+00 


-1.74E+00 


1.61E+01 


-2.19E+00 


-1.24E-01 


-3.45E+00 


1.9E-02 


cf 276 98 





-1.50E+01 


00E+00 

\J . V./V./ 1 — . \ \J\J 


-4.72E+00 


3.79E+01 


-4.15E+00 


-3.04E-02 


-9.12E+00 


3.8E-02 


cf 277 98 





-5.42E+00 


0.00E+00 


-3.49E+00 


2.64E+01 


-2.77E+00 


-1.12E-01 


-6.64E+00 


3.1E-02 


Example of Table 11 - 


— reverse (y, 


n)-rate fits on the basis of ETFSI mass-model predictions 






mother(A,Z) 


hu 


a rev 

"o 


a T 








a™ 


a T 




cf 269 98 





1.05E+01 


-3.46E+01 


-2.94E+00 


3.00E+01 


-4.69E+00 


1.40E-01 


-4.54E+00 




cf 270 98 





4.28E+01 


-5.82E+01 


-6.95E-01 


-3.73E+00 


1.60E-01 


-3.13E-01 


2.26E+00 




cf 271 98 





4.72E+00 


-4.43E+01 


-3.86E+00 


3.56E+01 


-5.02E+00 


1.90E-01 


-6.27E+00 




cf 272 98 





1.61E+01 


-5.62E+01 


-3.20E+00 


1.94E+01 


-2.75E+00 


-3.24E-03 


-2.99E+00 




cf 273 98 





8.14E+00 


-3.73E+01 


-4.68E+00 


4.14E+01 


-5.82E+00 


1.64E-01 


-7.89E+00 




cf 274 98 





2.49E+01 


-5.28E+01 


-2.67E+00 


1.32E+01 


-1.60E+00 


-1.64E-01 


-1.77E+00 




cf 275 98 





6.30E+00 


-3.57E+01 


-5.35E+00 


4.34E+01 


-4.94E+00 


2.04E-02 


-8.87E+00 




cf 276 98 





2.47E+01 


-5.01E+01 


-1.74E+00 


1.61E+01 


-2.19E+00 


-1.24E-01 


-1.95E+00 




cf 277 98 





1.00E+01 


-3.66E+01 


-4.72E+00 


3.79E+01 


-4.15E+00 


-3.04E-02 


-7.62E+00 




cf 278 98 





1.68E+01 


-5.02E+01 


-3.49E+00 


2.64E+01 


-2.77E+00 


-1.12E-01 


-5.14E+00 




Exampl 


i of Table 15 


— neutron-induced fission rate fits on the basis of ETFSI mass-model predictions. 




mother(A,Z) 


'fit 


ao 




a 2 








a 6 


Dev 


cf 268 98 





3.13E+01 


-6.04E-02 


8.24E+00 


-1.81E+01 


2.22E-01 


7.09E-02 


8.65E+00 


5.2E-03 


cf269 98 





3.38E+01 


0.00E+00 


2.22E+00 


-1.59E+01 


1.39E+00 


-7.64E-02 


4.52E+00 


4.5E-03 


cf 270 98 





9.75E+00 


3.15E-03 


-1.49E+00 


1.34E+01 


-1.77E+00 


1.59E-01 


-2.78E+00 


7.7E-03 


cf 271 98 





3.84E+01 


-5.73E-02 


8.81E+00 


-2.75E+01 


1.86E+00 


-9.37E-02 


1.04E+01 


1.5E-03 


cf 272 98 





-1.30E+01 


0.00E+00 


-4.24E+00 


4.22E+01 


-5.66E+00 


4.90E-01 


-9.00E+00 


1.7E-02 


cf 273 98 





3.87E+01 


-5.66E-02 


8.78E+00 


-2.78E+01 


1.89E+00 


-9.66E-02 


1.04E+01 


1.2E-03 


cf 274 98 


1 


-4.50E+01 


0.00E+00 


0.00E+00 


6.69E+01 


-8.82E+00 


6.52E-01 


-8.66E+00 


1.5E-03 


cf 274 98 





3.45E+01 


0.00E+00 


0.00E+00 


-3.97E+01 


3.13E+01 


-1.60E+01 


3.53E+00 


1.5E-03 


cf 275 98 





8.34E+00 


1.13E-02 


-2.48E+00 


1.50E+01 


-1.35E+00 


8.30E-02 


-3.67E+00 


1.1E-03 


cf 276 98 


1 


-3.74E+01 


0.00E+00 


0.00E+00 


5.39E+01 


-8.74E+00 


6.38E-01 


0.00E+00 


1.5E-02 


cf 276 98 





8.18E+00 


0.00E+00 


0.00E+00 


-2.88E+00 


6.77E+00 


-2.68E+00 


0.00E+00 


1.5E-02 


cf 277 98 


1 


8.74E+01 


0.00E+00 


0.00E+00 


-7.63E+01 


3.03E+00 


-7.28E-02 


3.10E+01 


1.2E-03 


cf 277 98 





1.87E+01 


0.00E+00 


0.00E+00 


-4.86E+00 


4.23E+00 


-1.66E+00 


4.15E-01 


1.2E-03 



Sec. 14.21 ). The value A' has to be multiplied by this ratio to derive 
the photodisintegration rate A: 

A y = A'-^- = A'x 1.18 = 5.9 x 10 9 s -1 . 

G274 Cf 

The values of the partition functions at Tc, are also given in the 
online Tables 3-6. Note that the procedure is always the same 
as described above, regardless of whether it is an exoergic or an 
endoergic reaction. 
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